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Preface 


This  work  on  eye-level  blood  flow  responses  to  +Gz  stress 
is  an  extension  of  the  work  performed  at  USAFSAM  in  1976  which 
modeled  eye-level  blood  pressure  versus  +Gz  stress'^.  At  that 
time  excellent  predictive  accuracy  was  obtained  in  relating 
mathematically  derived  blood  pressure  responses  to  actual  blood 
pressure  data  obtained  during  +Gg  stress  on  relaxed  humans.  At 
the  same  time  the  blood  pressure  data  was  being  generated  by 
placing  an  18-gauge  catheter  in  the  left  radial  artery , a 
Doppler  blood  flowmeter  was  placed  on  the  surface  of  the  skin 
over  the  temperal  artery.  Therefore,  data  for  the  previous 
work  on  blood  pressure  responses  and  that  for  this  report  come 
from  the  same  centrifuge  studies. 

Assistance  in  the  data  interpretation  was  gratefully 
received  from  K.  K.  Gillingham,  M.D.,  Ph.D.  of  USAFSAM  and 
S.  A.  Rositano,  Ph.D.  of  NASA/AMES.  The  generation  of  the 
analytic  transfer  function  from  the  ensemble-averaged  transfer 
function  was  accomplished  by  R.  C.  McNee,  M.S.  of  USAFSAM  on 
his  program  entitled,  "A  Fortran  Program  to  Fit  a Complex- 
Valued  Transfer  Function." 
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Materials  and  Methods 


Four  healthy  males,  all  members  of  the  USAFSAM  Accelera- 
tion Stress  Panel,  ranging  .in  age  from  24  to  26  volunteered 
for  the  studies  which  provided  the  data  analyzed  in  this 
report.  Doppler  ultrasonic  flowmeters  were  placed  over  the 
right  and  left  temperal  artery.  The  output  flow  velocities 
were  recorded  on  magnetic  tape. 

The  subjects  were  in  a relaxed  condition  seated  (13° 
setback  angle)  on  the  centrifuge.  Three  different  G stress 
profiles  were  taken:  1)  gradual  onset  runs  (GOR)  with  an 
acceleration  rate  of  0.067  G/sec.;  2)  rapid  onset  runs  (ROR) 
with  an  acceleration  rate  of  1 G/sec.  and  then  maintained  at 
a predetermined  level  for  either  15  or  60  seconds,  and  3)  a 
simulated  aerial  combat  maneuver  (SACM)  with  a total  period 
of  100  seconds  and  acceleration  ranging  from  0 to  1 G/sec. 
Profiles  were  chosen  to  challenge  the  subjects  visual  func- 
tioning without  causing  unconsciousness. 

Each  subject  was  told  to  be  in  a relaxed  state  and 
exposed  to  a GOR  to  the  visual  endpoint;  15  second  RORs  and 
at  least  one  60  second  ROR;  and  from  four  to  eleven  SACMs. 

At  least  one  60  second  ROR  and  one  SACM  elicited  a visual 
endpoint  on  the  subjects. 

Discrete,  finite  Fourier  Transforms  of  input  G stress 
and  output  (eye-level  blood  flow)  were  obtained  with  a 
Hewlett  Packard  5451A  Fourier  Analyzer.  To  prevent  spectral 
aliasing  and  also  convert  the  instantaneous  blood  flow  to 
mean  blood  flow,  both  input  and  output  signals  were  low-pass 
filtered  with  a 0.5  Hz.  48-db/octave  Butterworth  filter. 

Each  G stress  profile  contained  a 200  second  record  which 
was  sampled  at  .78  seconds.  The  spectra,  therefore,  had  a 
maximum  frequency  of  0.64  Hz  with  a resolution  of  5 mHz. 
These  sampling  parameters  were  used  for  the  resulting  reso- 
lution and  for  compatability  with  the  blood  pressure 
study  (1) . 


The  filtered  inputs  and  outputs  were  transformed  to  the 
frequency  domain  by  the  Fourier  Analyzer  and  the  individual 
transfer  functions  were  generated  by  dividing  the  output  eye- 
level  blood  flow  (0(f))  by  the  +G  stress  forcing  function 
(1(f)).  The  ensemble-average  of  the  G stress  runs  for  each  of 
the  four  individuals  and  the  ensemble-average  of  all  28  G 
stress  runs  was  calculated.  The  averaging  process  gave  equal 
weight  to  each  run.  The  rectangular,  polar,  and  impulse  re- 
ponse  of  the  ensemble-averaged  transfer  functions  of  each 
individual  and  the  total  ensemble-averaged  transfer  functions 
are  given  in  Appendix  A.  The  first  25  frequency  channels  (5mHz 
per  channel)  are  shown  for  the  transfer  functions.  The  first 
5G  channels  (.078  seconds  per  channel)  are  shown  for  the 
impulse  responses.  Only  the  first  25  frequency  channels  were 
used,  since  99.7%  of  the  +G  stress  forcing  energy  was  below 
.125  Hz  (25  channels) . 

An  analytic  transfer  function  was  then  generated  which 
most  closely  duplicates  the  frequency  profile  of  the  total 
ensemble-averaged  transfer  functions.  A Fortran  program  was 
used  to  fit  the  complex-valued  ensemble-averaged  transfer 
function  (2).  The  analytic  transfer  functions  sought  were: 

1)  triple  pole 

2)  double  zero  - double  pole 

3}  single  zero  - double  pole  - delay 

4)  single  zero  - triple  pole  - delay 

5)  double  pole  - delay  - 90°  rotation 

Parameter  estimates  were  obtained  for  each  analytic  term  by 
minimizing  the  weighted  sum  of  the  squared  error  for  both  real 
and  imaginary  terms  over  the  first  25  frequencies.  Marquardts * 
algorithm,  combining  the  Taylor  series  expansion  and  negative 
gradient  methods,  was  used  in  the  computations. 

Predictions  of  eye-level  blood  flow  responses  to  GOR, 

ROR,  and  SACM  G stresses,  based  on  the  actual  response,  the 
mean  actual  response  of  the  individual  and  the  response  gen- 
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erated  by  the  analytic  transfer  function,  with  the  least  stan- 
dard error  was  obtained.  Such  comparisons  gave  visual  evidence 
of  the  bidirectional  aspects  of  the  transfer  functions  and  the 
range  of  response  error  involved. 

Results 

The  mean  transfer  function  for  all  28  G stress  runs  is 
given  in  Appendix  A.  The  relatively  smooth,  well-defined,  con- 
tours of  most  of  the  empirically  determined  transfer  functions 
enabled  the  use  of  simple  analytical  estimates  of  a mathemati- 
cally derived  function.  The  various  equations  chosen  for  a 
possible  fit  with  the  standard  error  are  shown  in  Table  1.  A 
mathematical  derivation  relating  G stress  to  blood  pressure  and 
blood  flow  was  done  in  order  to  justify  the  analytic  functions. 
This  derivation  is  shown  in  Appendix  B. 

As  can  be  seen  from  Table  1 and  the  derivation  in  Appendix 
B,  the  analytic  function  which  best  fits  the  ensemble-averaged 
transfer  function  should  and  is  a gain,  delay,  differentiator, 
and  three  poles.  A comparison  of  this  analytic  transfer  func- 
tion and  the  one  generated  for  the  relationship  of  eye-level 
blood  pressure  versus  +G^  stress  was  also  accomplished. 

The  analytic  function  with  the  best  fit  for  eye-level 
blood  pressure  as  a function  of  +Gz  stress  was: 

H(s)  = -18.3  1 + 7.-I8s  +.  4--3-9.r^ 

1 + 3.99s  + 7. 93s^ 


The  analytic  function  with  the  best  fit  for  eye-level  blood  flow 

as  a function  of  +G  stress  was 

z 


H(s) 


-2.79s 

- ,936s  e 

1 + .741s  - 2.68s2  - 29.9s3 


This  would  mean  that  the  relationship: 

BF  _ delay  - differentiator 
BP  3 poles 


would  be  the  simplest  relationship  of  eye-level  blood  flow  and 
pressure.  This  relationship  also  is  logical  for  pressure  and 


Analytical  Function 


Standard  Error 


- .103 

2 5 

(1  - 1.80s  t 2.71s  - 15.2s  ) 


.0477 


- ,539s  (1  - ,0718s)  .0365 

1 + 3.04s  + 11.8s2 

-.115s 

- .537s  e .0364 

1 + 3.04s  + 11.7s2 


j(j  - 1.39W) 

- .119  e .0347 

1 + 2.09(ja)}  + 6 , 12  ( joj)  2 


-2.79s 

- ,936s  e 

1 + ,741s  - 2.68s  - 29.9s3 


,0295 


Table  1 


Analytical  Functions  Used  as  Estimates  of  the 
Transfer  Function  Relating  +Gz  Stress  to  Eye-Level  Blood  Flow 
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flow  in  an  artery  as  shown  in  Appendix  B.  The  differentiator 
and  delay  would  also  occur  since  different  arteries  were  used 
for  the  two  studies.  Thus,  the  differentiator  and,  especially, 
the  delay  could  be  due  to  a combination  of  the  theoretical 
derivation  and  data  collection  procedure. 

A mean  square  error  analysis  of  the  analytic  transfer 
function  (TFreg)  and  the  ensemble-averaged  transfer  function 
for  each  individual  was  also  accomplished.  This  analysis  is 
shown  in  Table  2.  This  analysis  was  done  on  a channel  (5mHz) 
basis  and  a total  square  error  basis.  The  analysis  showed  that 
individual  number  3 had  a drastically  different  transfer 
function  frequency  profile  from  the  other  three  volunteers. 
Therefore,  in  the  final  results,  the  analysis  of  the  blood  flow 
generated  by  +Gz  stress  profiles  for  volunteer  number  3 was 
neglected.  This  analysis  also  showed  aberrant  points  of  vol- 
unteer number  4 in  channel  number  5,  and  its  subsequent  second 
and  third  harmonics.  Since  these  frequencies  had  a abnormally 
large  error,  further  analysis  of  the  source  of  those  frequen- 
cies was  done.  It  was  noted  that  a transient  on  the  G stress 
curve  was  the  maximum  contributor  to  those  points.  Therefore, 
it  was  smoothed  to  compare  more  favorably  with  the  other 
ensemble-average  transfer  functions. 

Since  the  usefulness  of  the  analytic  transfer  function  is 
not  in  how  accurately  it  resembles  the  ensemble-averaged 
transfer  function,  but  in  how  well  it  predicts  actual  blood 
flow  responses.  We  compared  the  actual  response,  the  response 
predicted  by  the  individual  subjects'  mean  ensemble-averaged 
transfer  function,  and  the  result  predicted  by  the  analytic 
function.  These  curves  are  shown  in  Appendix  C. 

Data  interpretation  of  the  responses  was  done  visually 
and  by  the  HP5451A  Fourier  Analyzer  Unit.  The  two  main  points 
that  are  easily  visualized  are:  1)  the  erratic  actual  re- 
sponse of  subject  number  3,  and  2)  the  bidirectional  flow 
characteristic  of  both  the  analytic  transfer  function  and  the 
ensemble-averaged  transfer  function. 
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Due  to  the  erratic  behavior  of  the  flow  responses  of  sub- 
ject number  3,  his  actual  strip  chart  recording  data  was  again 
analyzed  visually.  It  was  noted  that  in  most  of  his  runs, 
blood  flow  and  blood  pressure  increased  before  the  decrease  in 

+G  stress.  This  could  be  explained  if  the  subject  did  not 
z 

remain  in  the  relaxed  state  for  the  entire  run.  Another  pos- 
sible cause  of  this  erratic  flow  could  be  the  state  of  the  art 
of  the  Doppler  flowmeter  at  that  time  (1974) . The  transducer 
and  signal  analysis  of  the  system  had  yet  to  be  optimized. 
These  erratic  responses  were  not  noted  in  the  previous  blood 
pressure  study,  since  they  did  not  affect  the  blood  pressure 
as  drastically.  This  could  be  due  to  increased  sensitivity 
of  the  blood  flow  data.  That  is,  percentage-wise,  the  blood 
flow  data  changed  more  than  the  blood  pressure  with  similar 
changes  in  -s-G  stress. 

The  bidirectional  flow  aspect  of  the  predicted  responses 
were  changed  to  unidirectional,  while  allowing  for  some  nega- 
tive flow,  by  inserting  a weighted  amplifier  on  the  output  of 
the  HP5451A  analyzer.  This  weighted  amplifier  is  shown  in 
Figure  1.  . u 


7 ‘ti  Tj 


out 


T*  4 ^ 

-C  kjl±.  tz 


Circuit  for  Changing  Bidirectional  Flow 
to  Weighted  Unidirectional  Flow 


The  dashed  line  on  +G  stress  runs  4506,  10936,  10530, 

z 

13100,  2300,  3906,  and  2730  in  Appendix  C are  generated  with 
the  weighted  amplifier  to  simulate  unidirectional  flow. 

The  time  error  analysis  was  done  on  the  HP5451A  Fourier 
Analyzer.  The  analysis  compared  the  time  differences  between 
the  actual  response  and  the  predicted  response  from  the  ana- 
lytic transfer  function  responding  to  the  28  +Gz  stress 
profiles.  The  maximum  and  minimum  points  on  the  responses, 
along  with  the  initial  point  of  zero  flow,  were  compared. 

The  results  of  this  analysis  are: 

Time  Differentials  for  Maximum  and  Minimum  Points 
Mean  = 1.5  seconds 
Maximum  = 4.5  seconds 

Time  Differentials  for  the  Onset  of  :ero  Flow 

Mean:  1.2  seconds 

Maximum:  5.0  seconds 

The  mean  time  differences  between  the  actual  response  and 
the  response  generated  by  the  analytic  transfer  function  are 
well  within  useful  and  anticipated  errors. 

Discussion 

Linear  systems  analysis  techniques  were  used  to  model 
nonlinear  systems.  This  limitation  was  deemed  minimal,  since 
the  ensemble-averaged  transfer  function  for  the  28  G stress 
profiles  correlated  well  with  the  delay-differentiator  three 
pole  analytic  transfer  function.  The  correlation  was  compa- 
rable to  that  of  the  blood  pressure  -G  stress  analytic 
transfer  function  (1) . 

The  blood  flow  data  was  generated  with  a Doppler  ultra- 
sonic system  which  cannot  yield  accurate  quantitative  data 
due  to  the  various  incident  phase  angles  encountered.  The 
results,  however,  are  modeled  for  percent  magnitude  changes, 
and,  therefore,  numerical  values  are  not  needed. 


-9- 


The  data  presented  in  this  report  has  been  theoretically 
compared  to  the  previous  work  on  blood  pressure  response  to  G 
stress  profiles.  There  is  no  other  work  which  could  be  used 
for  comparison  that  includes  blood  flow  data  as  a function  of 
G stress. 

The  time  analysis  of  the  data  shows  a mean  time  delay 
between  the  actual  response  and  the  analytically  generated 
response  of  1.5  seconds  for  maxima  and  minima  and  1.2  seconds 
for  the  onset  of  zero  flow.  It  has  been  shown  that  gray  out 
occurs  approximately  4 to  10  seconds  after  the  onset  of  zero 
flow  (6).  Thus,  the  mean  error  is  well  within  the  limits  for 
predicting  gray  out.  The  maximum  time  delay  of  5 seconds, 
between  actual  and  predicted  onset  of  zero  flow,  is  marginal. 
Thus,  prediction  capability  may  be  impaired  at  the  upper 
limits  of  error. 

Since  the  system  models  eye-level  blood  flow  of  relaxed 
humans,  it  can  only  be  a part  of  a real  world  model  relating 
G stress  to  eye-level  blood  flow.  A realistic  model  must 
incorporate  straining  maneuvers  and  the  use  of  anti-G  suits 
on  the  cardiovascular  system  dynamics.  Since  these  phenom- 
ena are  highly  nonlinear,  these  effects  should  be  modeled  . 
separately  and  then  added  to  the  basic  blood  flow  control 
model.  This  would  give  a complete  model  of  the  physiological 
effects  of  G stress  on  humans.  This  data  then  could  be  used 
effectively  to  predict  pilot  performance  during  various 
aerial  combat  maneuvers. 

Conclusion 

Data  describing  the  eye-level  blood  flow  response  to  G 
stress  in  relaxed  humans  has  been  presented  as  an  empirical  G 
to  blood  flow  transfer  function  and  as  a mathematically  well- 
defined  analytic  function.  The  transfer  functions  cover  the 
frequency  spectrum  from  5 to  200  mHZ  or  99.7%  of  the  forcing 
energy.  The  transfer  functions  exhibit  adequate  predictive 
ability  within  the  forcing  function  parameters  (1.0  -*■  4.5G; 

0 -*■  200  seconds)  . The  transfer  functions  were  modified  by 


an  operational  amplifier  circuit  to  enhance  the  unidirectional 
flow  with  weighted  reverse  flow.  The  delay-differentiator 
three  pole  function  most  closely  fits  the  empirical  transfer 
function  and  predicts  eye-level  blood  flow  with  usable  accuracy. 
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APPENDIX  A 


The  frequency  characteristics  of  the  various  transfer 
functions  used  are  given  in  rectangular  and  polar  coordinates. 
The  frequency  coordinates  are  given  in  the  left  column.  The 
frequency  value  is  equal  to  the  column  number  multiplied  by 
.005  HZ. 

The  "SF"  code  is  translated  below 
SF  ABC 
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APPENDIX  B 
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Theory  Guiding  Pressure-Flow  Relationship: 


The  steady  state  flow  is  described  by  Poiseuille  Law.  In 
the  transient  case,  it  is  essential  to  attempt  to  formulate  a 
precise  set  of  equations.  The  nonlinearites  have  to  be  consid- 
ered. Such  nonlinearites  arise  from  the  tapering  of  arteries, 
elastic  wall  of  blood  vessels  which  expand  due  to  pressure, 
etc.  (3).  Considering  these,  the  equations  for  pressure  flow 
based  on  Navier-Stoke  and  the  continuity  equations  are:  (4) 
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where 

S « area  which  depends  on  t,  p,  Z 
p = pressure 
t = time 

Z = axial  direction 

v = velocity 

ip  * term  for  seepage 

f = gravity  or  centrifugal  and  friction  forces 

It  has  been  shown  (3)  that,  with  friction,  etc.,  the  equa- 
tions can  be  written  as 


3 / c 2 , PS  . 

SZ  ( S v + T ) 


where 


S - 
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f » friction 
g « gravity 


TwD  , P 3S  3 ^ 

— + p 3z  = at  (Sv)  + pSg  {3) 
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D and  h are  diameter  and  wall  thickness  at  zero  pressure 
o o 

conditions. 

The  continuity  equation  is 

|f  ♦ ^ (vS)  =0  (4) 


If  flow  = 6 = vA,  then  it  can  be  shown  by  algebraic 
manipulations  that  the  Navier-Stoke  and  continuity  equations 
take  the  form: 
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These  equations  are  the  general  form  for  Navier-Stoke  and 
continuity  relationships.  A solution  of  these  nonlinear  partial 
differential  equations  must  be  done  by  approximation  techniques . (5) 
It  should  be  noted  that  pressure  and  flow  are  dependent 
not  only  on  g,  but  on  other  factors.  An  equation,  therefore, 
which  attempts  to  relate  p(t)  and  g(t)  or  0 (t)  and  g(t) 
directly  are  not  real  system  equations.  The  system  equations 
are  nonlinear  and  given  by  the  previously  nonlinear  equations. 

Algebraic  manipulations  of  the  Navier-Stoke  and 
continuity  equations  yields  the  following  differential  equations: 
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The  form  of  equation  which  relates  these  differential  equations 
is: 


3 ZP  , ! df  , , T3  _ /a  -7  6 IS  j a 17 
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36 


+ <J>4  Z6  } + {3^  + 32z  ||  + a3zg}  (8) 


where  A,  <j>f  y are  constants. 

The  above  equation  yields  the  general  form  of  the  blood 
pressure,  blood  flow,  and  gravitational  force  relationships. 
Note  that  pressure  and  flow  are  third  order  systems,  while 
gravitational  force  is  second  order.  This  shows  that  the 
transfer  functions  may  be  of  the  form 

„ -as 
Ke 

1 + as  +bs*"+"csJ 
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APPENDIX  C 


The  eye-level  blood  flow  responses  are  given  for  +G 

z 

stress  profiles  for  subjects  1,2,3  and  4.  The  four  curves 
are  arranged  in  the  following  order. 


FIRST 


+G„  stress 
z 


SECOND 


Actual  eye-level  blood  flow 
response  to  the  stress 


THIRD 


Response  generated  by  the 
individuals'  ensemble-averaged 
transfer  function  to  the  +G 
stress 


FOURTH 


Response  generated  by  the 
analytic  transfer  function  to 
the  +G  stress 


The  dashed  lines  on  the  responses  due  to  the  ensemble- 
averaged  transfer  functions  are  due  to  the  addition  of  the 
operational  amplifier  which  limits  reverse  flow. 
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